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0^ ■ We solve the problem of a magnetic impurity in a d^2_y2-w&Ye superconductor by a variational 

0^ ' method. A moment is found to exist in the superconducting state only if the Kondo-temperature 

in the normal state is larger than the maximum in the superconducting gap-function. If a moment 
exists in the superconducting state, and provided spin-orbit coupling is non-zero, it induces a time- 
' reversal breaking superconducting state locally around the inpurity which is a linear combination 

of the d^2_y2 and d^y states. The current pattern around the impurity in this state are evaluated. 
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INTRODUCTION 



o 

' Since magnetic impurities break time-reversal invariance, they tend to destroy superconductivityS. In superconduct- 
Oh' ing ground states paired in a finite angular momentum (or its appropriate generalisation in a lattice), non-magnetic 
^ , impurities_^e also pair-breaking since their potential, in general, does not transform in the same way as the order 
^ ^ • parameterB'El. 

, Recently a new phenomenon associated with magnetic impurities in a d^2_y2 condensate was proposed. It was 

• argued that the interaction between the magnetic impurity spin, and the orbital moment of the condensate can help 
I \ to stabilize a new time-reversal violating phase: adx2-y2+i(3dxya. The physical point is that such a state has a finite 

■ O ■ orbital moment around the impurity. Provided the spin-orbit coupling is finite, such a state interacts linearly with the 
S ■ magnetic moment. Therefore it is necessarily induced with the ratio P/a of order Eso/SEc where Ego is the spin-orbit 
coupling energy and SEc is the difference in condensation energy of the dx2_y2 and the dxy states. 

In this paper we examine this idea through a variational method introduced for the Kondo (and mixed- valence) prob.- 
lem in normal raetalsu. The variational method -foreshadowed the development of the non-crossing approximationQ, 
the 1/N methodQ and the Slave-Boson approachEl. For the ground state properties these methods produce the same 
results as the variational approach. First we ask the question: Under what condition can a moment exist in the 
superconducting state. This is the question of the disappearance of the Kondo-effect in a d-wave superconductor. 
This problem also has been examinedBEj. We hope , our simple approach, with results equivalet to those derived 
, earlier, makes the physical issues clearer. Because the density of states of quasi-particles goes to zero linearly in 
' energy, the logarithmic Kondo singularity in the scattering matrix of the moment with conduction electrons is ab- 
0^ , sent but, with a finite exchange coupling, the Kondo-effect and the disappearance of the magnetic moment is still 
possible. The condition for the moment to exist in the superconducting state is found to be Tr-/A ^ 1, where Tr- is 
the Kondo temperature in the normal state and A is the, nyximum in the superconducting gap-function. (This is a 
more stringent condition than in s-wave superconductors.Ej'tj). Second, we find the necessary condition for inducing 
the locally time-reversal breaking state. In agreement with earlier conclusions, a finite spin-orbit coupling is found 
. essential. With the variational wave-function, we also calculate the current distribution induced around the impurity. 

• We also briefly discuss the problem of several impurities to conclude that a global time-reversal breaking phase is 
Q , unlikely to result from this mechanism. 
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■ We consider the usual BCS Hamiltonian (Hbcs) and the Anderson Hamiltonian {Hj) to describe the superconductor 
and the impurity respectively: 

Hbcs - J2 ^k7L7k. + Eg (1) 

kcr 

with Ak = Y^e^ + , where Ck are the energy of conduction electrons with respect to the chemical potential, 

Ak = ^fefc'"^(c-kiCkT) is the superconducting order parameter, V^^F'^ is the BCS attractive interaction and c^^ 
creates an electron with momentum k and spin |. 



THE MODEL 
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7kT = -"kCkT - ■f^kcLkp 7-ki = WkC^T + "kC-ki (2) 
annihilate and create the quasiparticles in the superconducting state and 



The ground state is given by \G) — nk(^k + '^'kcl^^c^^i^^)\0) with an energy Eq- 
The impurity Hamihonian is: 



Hi ^ eo'^ d-l^doa + C/o^OT^oi + ^ Vk^o-Cka + H.c. (4) 



kCT 



djg. creates an electron with spin a in the impurity- orbital and no^ = dJ|C?o|. We take the limit Uq — > oo. 
With the Bogoliubov transformation, Eq. (H), the hybridization term (third term in Eq. (4) is: 

Hd,c = ^k[4T("k7kT + ^^kTlki) + dlii-Vkll^ + Mk7-ki)] + H.c (5) 

kCT 

where we have assumed Vk = V_k- 

We also include the spin-orbit interaction between the spin of the impurity (s) and the angular momentum of the 
conduction electrons (L): 

Hs^o ^9 dr — (6) 



with s = dg^(TQ,^do^ (cr are the usual Pauli matrices) and L(r) = cl^hcra- 

and L~ mix states which are even with odd under reflection about the x-y plane. As odd states have zero 
amplitude in the x-y plane and L~ are irrelevant in two-dimensions, which is the case we consider with Copper- 
Oxide metals in mind. Lz scatters a k state to a k' states with |k| = |k'|. 

The spin-orbit interaction can therefore be rewritten as: 



Hs-o = gsz L^^ cLck'CT = gsz Y ^CT^k'CT 

k,k'fT k,k'(T 

= gSz Y [(7kT7k'T - 7lki7-k'i)("kUk' + ^'k'Uk') + 7kT7lk'i(wk^'k' - "fk^k') + H.c 

k,k'CT 

Lz anticommutes with the time-reversal operator (T). Therefore L^^ = —L^ ^ = {L^ ^)* = iIm{L^ For planes 

waves ^ ~ sm{(p — ip'), where R is the radius of the sample. We absorb the coefficient in Lz defining a 

coupling g' . Lz is invariant under rotations around z and changes sign under reflection over planes which contains the 
z. If the problem has square symmetry in the x-y plane then we have a C^v group symmetry and Lz transforms as 
the one dimensional representation . The other possible representation of the same group are the one dimensional 
Aig{s), Big {d.j;2_y2) 8.11(1 B2g {dxy) , s.Tid the double representation E {p^, Py)- 
The model Hamiltonian considered in this paper is 

H = Hbcs + Hi + Hs-o (7) 



VARIATIONAL CALCULATIONS 



We consider a variational wavefunction which is a spin-doublet, reflecting a moment in the ground state, as well as 
a variational wavefunction which does not have a moment in the ground state. 
For the doublet we choose the wave-function: 

\D T) = (rfJt + E "^7.^1 + 4t E «kk'7^T^-k'i)|G) (8) 

k kk' 
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The first two terms are mixed by tlie hybridization, while the last term, which is a singlet pair of quasiparticles 
coupled to the local moment is mixed with the first term by the spin-orbit coupling under certain condi tions. The 
U = oo constraint is obeyed. 

The variational functions ak, Okk' are determined by the condition S{{D\H\D) — Ei:){S\S)) = 0, where the doublet 
energy {Ed) is referred to the energy of the BCS state. We get 

ED = eo-Y. "kVkM^ " f H L^'^'^iK^' - KK'Wk' (9) 

k kk' 

One should notice that ak and Ckk' have the symmetry of Vk and L^-Vk respectively. In a pure d^2_y2 superconductor 
flkk' has d^y symmetry (as T A2g ®^Big ^ ^B2g)- In the symmetric s-wave Okk' = because the form factor 
Ff — Wk^k' — VkUki vanish for |k| — |k'|. It is also worth noting that momenta k, fc'with different signs of the 
coherence factor Vk give the largest contributions. The complete expresion of ak and Okk' and Ejj are given in the 
appendix. 

Consider next the state without a moment in the ground state. Hbcs and Hj commutes with total spin (S't, 
impurity+conduction electron spins) while Hs-o does not. Then the eigenstates of H cannot be classified by the 
total spin St- As a consequence the usual singlet solution {St — 0) is replaced by a state with total spin projection 
z, Sz = 0, we will continue calling this state the singlet for simplicity. 

The simplest variational wave-function with Sz = has the form 

|5z = 0) = (l+^/?k4 + ^Afkk'7iTA'i+ E ^-^TUW-^"i)\G) (10) 

k kk' kk'k" 

•^k — "OT '-ki "Oi 'kT 

^k = 4T7lki + 4nU 

Again the last term can be non^Zjero only due to spin-orbit scattering. The first three terms naturally arise in a model 
with hybridization and U = odld. 

The complete expresion of the coefficients of Eq. (10) as well as the expression for the ground state energy £'5 are 
given in the appendix. 

In a normal metal, due to the Kondo-effect, the singlet always has lower energy. In the limit eg << — 

W 

eo^-2V'poH^—^) (11) 

where po in the density of states at the Fermi level and W is the half-width band. We obtain the usual Kondo 
temperature 

Tk = Ed-Es^W exp(-^) (12) 

2poJ 

with J = V^/\eo\. 

For an s-wave superconductor (Ak — Aq) in the same limit. Eg ^ eq + Aflr-sfhilc E^ eq. As a consequence 
the low energy state evolves from a singlet to the doublet when -eo/|V^| growsOO. In both cases the contribution of 
spin-orbit term is 0. 

The low energy limit of the density of states of a superconductor with nodes (points in 2-D, lines in 3-D) at the Fermi 
level is p{uj) = pquj/ Aq. The logarithmic singularity leading to the Kondo-effect is replaced by a term proportional 
to Lijln{uj) in the superconductor. An approximate analytic expresion of binding energy obtained from the equations 
( p^ , p5| ) in the appendix is 

^ -2V'p^[Ej1_^ ln(4^i_^) + ^ + ln( ^)] (13) 

'^"^ Ao ^Ej^-Ej+uoJ Ao ^uo, + {Ei^E'l)'^ ^ ' 

where we have assume a superconducting density of states linear in energy up to lJc Aq, and constant from there 
to W, and we have aproximated the self-energy terms in Eq. (p5|), — ri(Ak — Es) — r2(Ak — Es), by Ef) — eo. Using 
Eq. ([1^), it it possible to re-express Eq. ([l3| ) as 

Tpd rpd I rpd rpd I A 

^V^ln(^^V^)-ln(^). (14) 
Ao Aq Ik 
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We conclude from equation (|lj) that the doublet has lower energy if Aq > Tk- 

In Figure (1) we show Eo and Es as a function of Aq for a dx2-y2 order parameter Ak = AoCOs(2(/3), where (f 
is the polar angle, and an s-wave superconductor, Ak — Aq. We assume a constant non- interacting density of states 
Pq ~ 1/2W, with W the half-band-width. Other parameters are eo/W = —0.4, V/W — 0.28. It can be seen that while 
E]j remain almost constant and equivalent for s and d, Eg grows as A increases. This reflects the displacement of 
the low-energy excited states to higher energies and is faster for the gapped s-wave {Eg). As a consequence there is 
a crossing of levels at a given Ac. In the inset the effect of the spin-orbit coupling is shown. Both Ef, and Eg takes 
advantage of this term gaining a similar amount of energy and leaving A^ unaffected. 

Figure (2) shows A^ as a function of |eo|/M^ for different values of J. It can be see that A^ aproximately scales 
with Tk = W exp{—W/ J). The quotient AJ?/Tr- lies in the range 0.8 — 2.4 while Tk and Aq vary by several orders 
of magnitude. 

Changes in the macroscopic superconducting density of states p{lu ^ 0] will change the above results. A finite 
concentration of impurities in a d superconductor induces a finite — Op and helps to stabilized the singlet , while 
a complex macroscopic phase such as dx-^-y'-^ + idxy has the opposite effect. 

We can see how the doublet modifies the superconducting condensate: 

(Da\c^lc^^\Da)^ _ aj + Ek' «kk' ^ ^ 

{Da\Da) 1 + Ek "k + Ekk' ^kk' 

It can be seen that the presence of the impurity diminishes the (k f, — k J.) pairing; this reduction is of the order 1/N, 
(N is the number of sites), for one impurity but is measurable for a finite concentration of impurities. For Vk = V, the 
largest values of occur at the nodes of Ak . Therefore these states loose a greater relative weight. This tendency 
agrees with the extended gapless region located in the gap function centered around the nodes, found in a momentum- 
dependent scattering for non-magnetic impurities in unconventional superconductorso. The singlet has a similar effect 
on the correlation function (c_kiCk|) ( a| + J2k' "^kk' should be replace by 2/3^ + M^^^, + X^k'k" ^kk'k" Eq.(13) 
). In this case a bigger distortion is expected (in general f3k is greater than ak). 
The spin-Orbit interaction introduces a new correlation: 

/n m \i^D /n I in \ i9<yIm{L^^' ){uk'Vi^' - VkUk) 

{D(j\D(j)F^^, = (Dcr|c_kiCk'T|£''7) = 7T — --; — (16) 

(Ak + Ak' + od) 

If the symmetry of the pure superconductor is dx2_y2, F^^, ^ id^y. We can calculate the induced order parameter in 
real space: 

A.fR - (-P'^lcR-riCR+rTlJ'g) _ ^ V- ,R(k-k') ir(k+k')pD f.^) 

{Da\Da) ^ 

Note that as k 7^ k', A^ depends on R. Taking the impurity site (R ~ 0) as the center of the symmetry this 
state has a dxy symmetry. Centered at a different site R, Ai(R, r) does not have a definite symmetry in the relative 
coordinate r. 

This state produces a complex pattern of spontaneous currents around the impurity which can be seen evaluating 
the current density operator: 

j{R){Da\Da) = {Da\j{R)\D<7) ^ dc ^ e'^C^-'^') (k + k')(uki;k' - w)a;;k' (18) 

kk' 

where the constant c = eh/mVol. 

Expanding the form factor (uki'k' ~^^kWk' )^ of Eq. ( |l8|) in the first two spherical harmonic we can get an approximate 
analytical form for this current, 

where J„(i?) is the usual Bessel function of order n. 

In Fig. (3) we show the current density j(R). The magnitude and direction of the arrows represent the current at 
that point. This pattern can be understood as the sum of a net current around the impurity and eight secondary 
currents. Four of them are around the points Rkp = (±4, 0) and Hkp = (0, ±4), and have the same sense as the net 
current, and the other are around the points Rkp = (±4, ±4), with the oppossite sense (see Fig. (4)). The sense of the 
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net current depends on the impurity spin projection Sz- Fig. (5) shows the net current flowing (J(i?) = ^ ^j(R)d(/5) 
as a function of the distance to the impurity. As before we can get the approximate analytical results, 

^(«)c< 4^+4^1™ + ,.^ (20) 

The sign of the current varies with distance but the total current around the impurity is not zero. We have neglected 
any diamagnetic current produced by the magnetic field induced by the spontaneous current. 

For the singlet the two momenta correlation gives F^^^, = M^k' (wk'itk — v\t_v\^i), but for Uk real (real macroscopic su- 
perconducting state), this state preserves the time- reversal symmetry {T\S) = \S)) and does not have any spontaneous 
current. 



CONCLUSIONS 



In summary we have shown, in agreement with previous results, that a magnetic impurity in a dx2_y2 superconductor 
has a transition from a Kondo singlet to an unscreened doublet as the superconductivity is turned on. The coupling 
between the spin of the impurity and the orbital momenta of the electrons induces a complex secondary component 
of the superconducting order parameter around the impurity in the doublet. 

What are the observable consequences of the complex order parameter around the magnetic impurity for the case 
of dilute impurity concentrations? A complex order parameter by itself leads to a finite gap in the excitation spectra 
of the supeconductor. But the potential scattering due to the impurity, not considered in this paper, if it is not sitting 
in a center of symmetry of the crystal, produces a finite density of states at low energies and spoils this effect. The 
principal effect of such a potential scattering in the current around the impurity is that it decays in a length of the 
order of magnitude of the mean-free path, rather than decay as a power law obteined in this paper. 

Global time reversal breaking due to the current around the magnetic impurities is possible only if they are aligned 
ferromagnetically. This is because the direction of the current depends on the direction of the magnetic moment. 
Ferromagnetic aligment is unlikely in general. A spin glass state, with an associated glassy pattern of currents, is 
much more likely. If a ferromagnetic aligment of the impurity spins is achieved, the nature of superconductivity is 
strongly affected-a|nd may be even destroyed. The physical properties in such a case have been discussed for s-wave 
superconductorstSIU 
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APPENDIX A 

The coefficient of the doublet (Eq. ||) are given by: 

■ -9 r /rkk'N/ ^ , VkUkVk'Wk' , ^k^k V^T^ /01^ 

flkk' = ~^77, Im{Lz )(ukWk' - WkUk'j + tt + tt 2^ 14"Wk"ak"k (21) 

^'kk' l^k - -C/D jtkk' l^k - jfkk' ^ 

'kk' = eo + Ak + Ak' — Ed 

^ -^kMk ^ g V^. v* Im{L^^') ~ '"^'^^'^ (22) 

Ak — Ed — Fak 2(Ak — Ed — F^k) ^ ^ luw 

"■^ eo + Ak + Ak' - Ed 

k' 

The second terms in Eqs. ^ and (ak and Okk') only gives a contribution to the energy if Fy F„ F^^ = Fy, 
where Fy, Vy and F^^ are the respective irreducible representations. In a Cay group this is only possible if Fy — E 
{Px, Py)- This means that the excited (second and third terms in Eq. pT] ) connected with via H, are not connected 
to each other. 

Considering the most relevant terms the doublet energy Ed is given by 
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ED-eo = -Sd = -Ti{-Ed) - r2{-ED) 



(23) 



with: 



ri(x) = 5] 

k' 



4 ^ eo + Ak' + Ak + X 



APPENDIX B 

For the singlet (Eq. the ground state energy is calculated to be: 

Es^eo-5s^2Y,l3^V^< (24) 

k 

where 

" Ak + eo - - ri(Ak - Es) - r2(Ak - Es) + ^^^^ 



The coefficients are given by 



,cor _ Ek' /3k' (^L;^. - ^Lk' + «5/2(i.(uu* + t'^;*)kk') 



" Xk + eo-Es-T^{Xi,~Es)-T2{Xk-Es) ^ ' 

(1) _ Vk^kVk'Mk' 

" Ak + Ak' - Es 



.(2) _ 5^ V 

^kk' ~ A 2^ 



(2) _ 5^ L^(ut> - t;M)k'k"^z(™ - t'M) 



kk" 



4 t;;' Ak + Ak' + Ak" +eo-Es 

-(Vk^k/3k' +Vk'Uk'/3k) 

= Ak + Ak^-i?. 

&kk'k" = „.g /3k/m(L^)(w - i;u)k'k" |,2g^ 
2 Ak + Ak' + Ak" + eo — Es 

^kk'' ^kk' ^'^^ ^'^^ ^^^^ term in Eq. ^ contribute only in very special cases (i.e. ^kk' gi'^^s a contribution only 
if Fy ® Fy = F„ (Fy = F^, = Aig{s) or Fy = E(px,Py)). They do not contribute for the d-wave superconductor 
consider in the main text. 
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CAPTIONS 

Figure (1) - Doublet energy, = Ef^ (full line), and singlet energies Eg (dotted line) and Eg (dashed line) as 
a function of Aq, for s-wave and d^^_y2 superconductors. A* and A^ denote the values of the superconducting order 

parameter whore the levels crosses. Other parameters are e^/W = —0.4, V/W = 0.28.. A,.. The inset shows Ef) and 
Eg with the inclusion of an spin-orbit interaction, with a constant coupling g' /W = .1, (dashed lines) compare with 
the g' = case (full lines), near the crossing A^. 

Figure (2) Critical values Af./TK{J) as a fimction of |eo|/W^ for different values of J. 

Figure (3) Current density j(R) around the impurity. The magnitude and direction of the arrows represent the 
current at that point. The sense of the currents depends on the impurity spin projection Sz- 

Figure (4) Current density j(R) around the impurity after substracting the net current (Fig. 5). The magnitude 
and direction of the arrows represent the current at that point. 

Figure (5) Net current J{R) = l/27r^ j(R)di^ as a function of the distance to the impurity. 
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